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Abstract 

We seek possible statistical consequences of the way a forcing term is added 
to the Navier-Stokes equations in the Direct Numerical Simulation (DNS) of 
incompressible channel flow. Simulations driven by constant flow rate, constant 
pressure gradient and constant power input are used to build large databases, 
and in particular to store the complete temporal trace of the wall-shear stress for 
later analysis. As these approaches correspond to different dynamical systems, it 
can in principle be envisaged that these differences are reflect by certain statistics 
of the turbulent flow field. The instantaneous realizations of the flow in the 
various simulations are obviously different, but, as expected, the usual one-point, 
one-time statistics do not show any appreciable difference. However, the PDF 
for the fluctuations of the streamwise component of wall friction reveals that 
the simulation with constant flow rate presents lower probabilities for extreme 
events of large positive friction. The low probability value of such events explains 
their negligible contribution to the commonly computed statistics; however, 
the very existence of a difference in the PDF demonstrates that the forcing 
term is not entirely uninfluential. Other statistics for wall-based quantities (the 
two components of friction and pressure) are examined; in particular spatio- 
temporal autocorrelations show small differences at large temporal separations, 
where unfortunately the residual statistical uncertainty is still of the same order 
of the observed difference. Hence we suggest that the specific choice of the 
forcing term does not produce important statistical consequences, unless one is 
interested in the strongest events of high wall friction, that are underestimated 
by a simulation run at constant flow rate. 
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1. Introduction 


A flow through a straight duct always requires an external force to overcome 
the friction losses at the wall. In a pipeline, this force is typically generated by a 
pump; in a drain pipe, it is provided by gravity. Similarly, numerical simulations 
of duct flows also require a modelling step where this external force is defined. If 
we restrict our attention to the direct numerical simulation (DNS) of turbulent 
channel flows, where the incompressible Navier-Stokes equations are considered 
a more than satisfactory mathematical model, the most popular choice for the 
driving force is implicit in the prescription of a Constant Flow Rate (CFR): the 
time-dependent force is calculated by adjusting the uniform (in space) pressure 
gradient required to keep the flow rate constant. Alternatively, a Constant 
(in time) Pressure Gradient (CPC) can be prescribed, thus explicitly setting 
the driving force, which becomes constant and uniform. For a flow through a 
straight duct where friction at the wall is the only source of energy dissipation, 
the constant value of the pressure gradient is proportional to the mean value 
of wall friction, and the simulation a posteriori yields as a result the mean flow 
rate required to generate the prescribed frictional losses. As a third alternative 
besides CFR and CPC, we have recently proposed a Constant Power Input 
(CPI) strategy 0, where the product of flow rate and pressure gradient is kept 
constant at every time step, while either one of them can instantaneously vary. 

For laminar flows these three simulation strategies must yield identical re¬ 
sults, since flow rate and pressure drop are uniquely related. An impact of the 
forcing strategy can only be seen when the effect of changes in the flow geometry 
is studied, and the flow field with and without this change is to be compared. 
Additional frictional losses will lead to an increased pressure gradient for CFR or 
a decreased flow rate for CPC. The resulting difference can nicely be visualized 
in the money-versus-time plane [ij . 

On the other hand, in turbulent flows the different forcing strategies lead 
to different time histories of the flow quantities. Such differences are generally 
considered to be without statistical consequences. It is generally known that the 
various forcing strategies are not equivalent on practical grounds; for example it 
is known [l^ that, when the value of the Reynolds number is suddenly changed 
in a simulation, CFR does present shorter transients and is faster to adapt to the 
new state of statistical equilibrium. This is a practical advantage that represents 
one of the reasons why the CFR approach is so far the most popular. However, 
the general consensus is that, once the statistical equilibrium is reached, the 
flow statistics do not depend upon the particular choice of the forcing term. 
In other words, as long as the same turbulent channel flow is considered in all 
cases, its statistics are independent from the forcing term. 

Nevertheless, in principle a difference exists between the forcing strategies 
in turbulent flows. The CFR strategy results in a temporally fluctuating pres¬ 
sure gradient, while CPG produces a flow rate that fluctuates in time, and 
in GPI both quantities fluctuate. These fluctuations depend on the computa¬ 
tional domain size and become smaller for larger domains Q, but are always 
present in a domain of finite size. Owing to this difference, it can be envis- 
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aged that differences in the flow dynamics might exist, which have so far not 
been identified in the typically considered statistical quantities. In the present 
contribution, we will specifically tackle this fundamental question, by compar¬ 
ing various statistical quantities computed from three purposely built databases 
for turbulent channel flow computed with CFR, CPG and CPI at nominally 
the same Reynolds number. The underlying data sets are very large to ensure 
that errors in the statistical estimates are as small as possible. The analysis is 
mainly focused on the wall-shear stress, which is a quantity that can be seen as 
a footprint of the entire turbulent flow field above the wall. 


2. Forcing strategies 


We consider a fully developed turbulent channel flow of an incompressible 
Newtonian fluid. The channel flow is set up in such a way that the stream- 
wise coordinate corresponds to Xi or x, the wall-normal one to X 2 or y and the 
spanwise one to X 3 or z. Throughout the entire paper dimensional quantities 
are marked with an asterisk, while unmarked symbols correspond to dimen¬ 
sionless quantities. A -|- superscript indicates quantities made dimensionless in 
inner units; otherwise, the length and velocity scales are tose employed in the 
definition of the Reynolds number. 

In order to drive the fluid through the channel, the Navier-Stokes governing 
equations have to be supplemented by a forcing term, i.e. a force per unit 
volume, so that they read: 


9uj 

dxj 

dui dui 

dt dxj 
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( 1 ) 
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dxi Re dxjdxj 
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where repeated indices imply summation, Ui represents the velocity component 
in the *-th direction of the Cartesian reference system, and Xi is the correspond¬ 
ing spatial coordinate. The Reynolds number Re is formed by the channel 
half-height S* and the characteristic velocity scale for each forcing strategy, i.e., 
the bulk mean velocity for CFR, the friction velocity u* for CPC or the 
power-based velocity for CPI. As explained in [^, is the bulk velocity 
obtained in a laminar flow with the imposed power input. The additional forc¬ 
ing terms in the three components of the momentum equation are denoted by 
A,. 

The present paper seeks for possible statistical consequences of selecting the 
streamwise component Fi of the body force, the definition of which in depen¬ 
dence of the three forcing strategies (CPC, CFR and CPI) is described below. 
In respect to the other two spatial directions it is obvious that F 2 = 0, as no flow 
rate or pressure difference can exist between the two impermeable walls. Very 
often it is also considered that F 3 = 0. Strictly speaking, the two homogeneous 
directions are equivalent, so the question about the possible relevance of the 
forcing term should be asked in regard to the spanwise direction, too. However, 
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the constant value (of the spanwise flow rate or the spanwise pressure gradient) 
enforced in the simulation is zero, hence possible effects, if any, are thought to 
be of a lesser entity. In this paper we restrict ourselves to considering Fi , which 
is labeled F from here on. 

2.1. Constant Pressure Gradient (CFG) 

The simplest forcing strategy is to set F to a constant value. The forcing 
term can be regarded as a spatially averaged pressure gradient: 



( 3 ) 


where (•) indicates the spatial average in the homogeneous directions x and z. 
Therefore, we refer to this strategy as CFG. Since the pressure gradient must 
balance the wall friction in a fully developed flow, the following force balance is 
satisfied: 



( 4 ) 


where is the mean value of the wall friction. By using the friction velocity 
K = where p* is the density of fluid and d* the channel half height, 

the dimensionless forcing term becomes F = 1. The resultant Reynolds number 
based on u* in Eq. ([2]) is referred to as the friction Reynolds number Rer, which 
is given by Rer = u*6* Jv*, where v* is the kinematic viscosity of fluid. 

2.2. Constant Flow Rate (CFR) 

In the CFR strategy, which is often employed owing of the practical advan¬ 
tage mentioned above, the forcing term F is not constant in time anymore, but 
changes at every time step in such a way that the resultant flow rate is equal 
to a prescribed flow rate. The flow rate per channel cross section is the bulk 
velocity and the corresponding Reynolds number for CFR is thus the bulk 
Reynolds number Reb = U^S*/v*. 

Although the CFR condition is commonly used, its detailed numerical im¬ 
plementation is quite often not fully described in the existing literature, and 
more than one variant exist which differ in the specific algorithm that keeps the 
flow rate constant. Sometimes the CFR strategy is detailed when the numerical 
method has to deal with additional complications: for example in Q CFR is 
explained in the context of incompressible and compressible simulations, while 
in 0 CFR is described for the simulation of an incompressible flow in a channel 
with riblets on the solid walls. 

Typically, the Navier-Stokes equations m is discretized as follows: 



where Sn is the Kronecker delta, a variable at a time step of n is denoted by the 
superscript n, and A/) represents all non-linear terms in the z-component of the 
momentum equation. As the full velocity held up to time step n is available, the 
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non-linear term TVj ^ is often calculated explicitly by using available informa¬ 
tion at the previous time steps, i.e. with explicit methods like Adams-Bashforth 
or Runge-Kutta. As for the viscous term, an implicit scheme such as the Crank- 
Nicolsoir irrethod is often used in order to mitigate the requirement for a small At 
to avoid numerical instability, which becomes crucial especially in the near-wall 
region, where fine grid resolution is generally employed. The pressure at time 
step n -I-1 is often calculated by the fractional step method. Otherwise, Eq. m 
may be solved by a wall-normal velocity - wall-normal vorticity formulation, 
where the pressure term is eliminated. In a simple flow geometry as the one 
considered here, when the boundary conditions are periodic, this formulation 
presents considerable advantages. It is important to notice, however, that the 
following discussion can be applied without loss of generality regardless of the 
numerical scheme used to solve Eq. discretization is introduced here only 
for establishing notation. 

In the streamwise direction the forcing term that is required to keep 

the flow rate constant has to be determined. However, the flow rate at time step 
n -|- 1 is only known as a result of the computation. Therefore, the flow field at 
time step n -I- I is assumed to be the superposition of two components, namely 


u 


n-\-l 

1 




- 72+1 


( 6 ) 


Here, represents the solution for the streamwise velocity component at 

time step n + 1 without application of the forcing term, i.e. is obtained with 
= 0. Substituting Eq. ([ 6 ]) into the first component of Eq. m results in 
the following equation for the second term, 


- 71+1 


At 



I 

2Re dy'^ J 


(7) 


Here, we use the fact that u is a function of y only, since the forcing term does 
not depend on x and z according to its definition ([3]). The analytical solution 
of Eq. ([7]) can be easily obtained by taking into account the no-slip condition 
at a wall {y = 0 ) as: 


~n+i ^ _ exp(-Ay)} , (8) 

where A = +2i?e/At. Note that this solution is essentially uniform away from 
the wall and only changes rapidly in the vicinity of the wall due to the viscosity. 
The increment of the bulk mean velocity due to is obtained by integrating 
Eq. ([5]) from the wall to the channel center: 

AH,"+i = = AtF"+i (^1 - . (9) 

The numerical procedure for CFR can thus be summarized as follows. First, 
is obtained by any numerical scheme with = 0. The resultant flow 
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rate, i.e. the wall-normal integral of the profile, now differs from the pre¬ 
scribed value. Hence, the flow rate is compensated by adding a corrective profile 
proportional to given by Eq. ([5]). The amplitude of is determined 
so as to satisfy the constant flow rate condition exactly based on Eq. (0, which 
can be solved only once if the DNS integrates the equations of motion with a 
constant time step, whereas a condition based on the CEL number implies a 
variable time step and the need for solving Eq. whenever the value of the 
time step is changed. 


2.3. Constant Power Input (CPI) 

Recently, the present authors have proposed Q a third simulation strategy 
that is alternative to the CPG and CFR conditions explained above. In this 
strategy, it is the power input to the flow system that is kept constant in time, 
instead of the pressure gradient or the flow rate. The corresponding character¬ 
istic velocity C/j^ is one that is based on the power input P* to the flow system 
per unit wetted area: 


and is given by: 



( 11 ) 


where /i* is the dynamic viscosity of fluid. This quantity corresponds to the 
bulk velocity achieved by a laminar flow under a given power input P*. In 
the case of turbulent flows, the resultant bulk mean velocity is smaller than 
C/n owing to the additional momentum loss induced by turbulent stresses. The 
corresponding power-based Reynolds number is defined as Ren = 11(^6*/v*. 

Under the CPI condition, the product — is kept constant according 

to Eq. (unj. This is expressed in dimensionless form as: 


Pp — 


(^\u = = — 

dx / ^ P*U))^ Ren 


( 12 ) 


so that in case of CPI the flow is driven by the fixed power based Reynolds 
number Ren- 

The implementation of the CPI condition to DNS is not trivial, as the power 
input at the next time step Pp^^ depends on the flow rate which is 

determined as a result of the computation. An expression for the corresponding 
forcing term can be obtained using the relation R = ~ such that (as 

introduced in Eq. ([5|)) can be approximated with first-order accuracy in time 
as: 


p 

pn+l _ fjl 

UII' 


(13) 


The increase in the order of accuracy is straightforward, and it has been con¬ 
firmed that further increase of the accuracy does not affect the flow statistics 
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forcing strategy nominal Re measured Rct 
CFR Reb = 3173 Rcr = 199.51 

CPG Rer = 200 Rcr = 199.97 

CPI i?en = 6500 i?e^ = 199.71 

Table 1: Values of Re for the three simulations. The nominal Re for each simulation is shown 
together with the value of Rct measured a posteriori. 


for the small time steps typically used in DNS of turbulent channel flow. The 
detailed numerical implementation of the CPI condition is described in [^. 

3. Results 

Three databases have been built by carrying out three DNS of a turbulent 
channel flow based on the three different forcing strategies described above. The 
DNS code is the one described in [l^ which employs mixed spatial discretization: 
Fourier series in wall-parallel directions and compact hnite differences in wall- 
normal direction. The flow is in all cases at a nominal Rcr = 200, and the very 
same discretization (in space and time) is used. The computational domain 
is 47r5* X 25* x 27:5* in the streamwise, wall-normal and spanwise directions, 
and matches that employed in the seminal DNS study by Kim et al The 
number of Fourier modes and grid points is = 256 and Ny = 128, 

with a spatial resolution of Ax~^ = 9.6, A?/+ = 0.8 — 4.9 and Az+ = 4.8. The 
spatial resolution in the homogeneous directions is further increased thanks to 
the complete removal of aliasing error via temporary increasing the number of 
Fourier coefficients by at least 50%. The time step is fixed at At+ = 0.2, and 
the calculations last for 150,000 viscous time units. 750 independent flow fields 
for each case are stored for computing one-time statistics. When computing 
statistics, the two halves of the channel are averaged together to double the size 
of the statistical sample. The main focus of the analysis is wall friction, the two 
components of which are written to file for further analysis every one viscous 
time unit. The complete time history of wall friction thus occupies 300 CB of 
disk space for each case. 

Table [T] lists, for each of the three cases, the numerical value of the specific 
Reynolds number that was given as input to the simulation, together with the 
value of Rcr that is measured a posteriori by time averaging over the full sta¬ 
tistical sample. It can be seen that all the cases achieve the target Rer = 200 
within an extremely small tolerance. In the CPC case Rer is specified directly as 
the input value, so that the tiny difference between the target and the measured 
value is simply ascribed to the time-averaging process. The non-CPC cases re¬ 
quired a little trial-and-error to identify suitable values of Reb and i?en, and 
this explains the relatively larger difference between the nominal and the actual 
value of Rcr- However, the present differences in Rcr are so small (0.25% in the 
worst case) that the simulations can safely be considered at the same nominal 
value of Re. The non-dimensionalization in inner units of the results presented 
in the following is based on the actual value of Ur obtained for each case. 
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Figure 1: Time histories of the bulk velocity Uh (dashed lines) and the forcing term F (contin¬ 
uous lines) with CFR (red), CPG (black) and CPI (green). Time-varying values are plotted 
over a relatively short time window (only 1/15 of the total length of the simulations) and are 
normalized by their long-term averaged values Ub and F. The red dashed line is not visible 
as it is constant and completely overlaps the black continuous one. 


Differences between CFR, CPG and CPI are visible when time traces of 
bulk velocity and space-mean streamwise pressure gradient (equivalent to the 
forcing term F) are compared. In general, the amplitude of the fluctuations 
depends upon the domain size and increases with smaller domain sizes |,i, 
being rather small but not negligible for the present case, where the domain 
is rather large as in [^. The forcing term itself fluctuates with CFR, the flow 
rate fluctuates with CPG, and both vary in time when a constant power input 
is prescribed in CPI. Figure [T] shows time traces of the bulk velocity and the 
streamwise pressure gradient for all three cases for an arbitrarily selected and 
relatively short time window (only 1/15 of the total length of the simulation). 
The quantities are normalized with their long-term average. It can be observed 
that fluctuations of the pressure gradient for CFR are characterized by larger 
relative amplitudes and higher frequencies than fluctuations of the bulk velocity 
for CPG. Fluctuations of the flow rate for CPI are qualitatively similar to the 
ones observed for CPG. The fact that the product of flow rate and pressure gra¬ 
dient is constant for CPI can be easily appreciated visually by noticing that the 
time histories for both quantities, i.e. the two green curves, are mirror images 
along the horizontal axis. We also note that the temporal fluctuations of the 
pressure gradient and the bulk velocity in CFR and CPG, or both in CPI should 
all vanish if the computational box is infinitely large. In practice, however, the 
computational domain size is always finite, and therefore the different features 
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Figure 2: Wall-normal profiles (in inner scaling) of the longitudinal mean velocity (left) and of 
the r.m.s. value of the fluctuations of the velocity components (right), for the three simulation 
strategies. In the right figure, the continuous line corresponds to the streamwise component, 
the dahsed line to the wall-normal one, and the dash-dotted line to the spanwise component. 
No difference can be appreciated. 


of the temporal oscillation in the forcing term may affect the flow dynamics and 
the resultant flow statistics. 

In contrast to time traces, which by themselves bear no statistical signifi¬ 
cance, the main first- and second-order single-point and single-time statistics, 
like the profiles of the longitudinal mean velocity or of the r.m.s. amplitude of 
velocity fluctuations, are expected not to present discernible differences among 
the three databases, in accordance to almost 30 years of practice in the DNS 
field. This expectation is fully confirmed. An example is shown in figure [21 
where the three different curves for CFR, CPG and CPI are indeed undistin- 
guishable from each other, as they overlap almost perfectly. The figure plots 
the mean streamwise velocity profile in linear scaling as well as the profiles of 
the r.m.s. values of velocity fluctuations. Similar results hold for other such 
statistics, among which we only show in figure [3] one-dimensional spectral den¬ 
sity functions for the velocity components in the first inner grid point, with the 
sole aim of showing their smoothness, made possible by the sheer size of the 
databases. 

For further analysis, we focus in particular on the statistics of streamwise wall 
friction, for which we have stored full time histories and not only uncorrelated 
flow fields. As indicated by the Reynolds numbers reported in Tabled! all three 
data sets have an almost equivalent friction Reynolds number of Rer = 200 
which indicates a nearly identical mean value of wall friction. This is confirmed 
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Figure 3: One-dimensional spectral density functions (inner units) for the velocity components 
at the first inner grid point, as a function of the streamwise wavenumber kx (left) and of the 
spanwise wavenumber kz (right). Different lines correspond to different velocity components 
while the forcing strategy is coded with different colors as indicated in figure [2] The red and 
black lines cannot be seen due to the almost perfect overlap of the data. 


forcing strategy r.m.s. third mom. fourth mom. 

CFR 0.3719 0.9405 4.2924 

CPG 0.3720 0.9380 4.2895 

CPI 0.3720 0.9400 4.2945 

Table 2: Statistical moments for the streamwise component of wall friction. The second 
moment is expressed as the root-mean-squared value of the fluctuations divided by the mean 
value. Third and fourth moments are by definition dimensionless. At Rcr = 180 values of 
0.367, 0.930 and 4.219 are calculated in Q- 


10 










Figure 4: Probability density function of the instantaneous, local fluctuations of the stream- 
wise component of wall friction around their time- and space-mean value. Abscissa is in units 
of standard deviation. Shown are the present results (left) and results taken from [3 at the 
slightly lower Re-r = 180 (right). For the present data negative local friction corresponds 
to O' < —2.79, marked by the vertical dashed line. No difference between the three forcing 
strategies can be appreciated. 


by the values shown in Table [21 where additionally the r.m.s. value of the 
fluctuations expressed in percentage of the mean is given, together with the 
dimensionless skewness and flatness factors. All of them agree very well among 
themselves as well as with available published data Q- 

It has recently been recognized that rare events in the near-wall flow dy¬ 
namics may lead to localized backflow, hence to negative values of streamwise 
friction. Lenaers et al Q made this observation in turbulent channel flow, but 
similar results are also observed in the zero pressure gradient turbulent bound¬ 
ary layer [1^ . By looking at the probability density function for the wall 

friction, it was determined in [7| that the probability of local backflow is around 
0.01% at Rer = 180, and has a tendency to increase with Re. Here we compute 
the same PDF and confirm that result; actually a slightly higher probability of 
about 0.02% for backflow is observed, which sounds reasonable given the slightly 
larger value of Re considered here {Rcr = 200 vs Rer = 180). OveraT the PDF 
for the three forcing types look very similar to the one reported by [7| , without 
revealing any difference related to the forcing strategy. 

A different scenario emerges when the same PDF is plotted by taking advan¬ 
tage of the large size of the database, that comprises several billions datapoints. 
Such a database allows extreme events to be tracked, and the PDF plot can be 
extended to contemplate much lower probability values. This full plot is shown 
in figure m where the PDF is observed to extend for almost nine decades. Much 
wider tails are evident now, with the tails being remarkably exponential (straight 
lines in this plot), especially for cr > 0. The PDF are obviously quite irregular 
in the extreme part of their tails, owing to the small number of events belong¬ 
ing there. Nevertheless, the forcing term, while apparently not affecting the 
extreme events of negative friction, does indeed appear to affect the strongest 
wall shear stress events with positive friction. In particular, CFR presents less 
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Figure 5: Probability density function of the instantaneous, local fluctuations of the stream- 
wise component of wall friction around their time- and space-mean value. The abscissa is in 
units of standard deviation. Negative local friction is found for cr < —2.79, marked by the 
vertical dashed line. The two oblique straight lines emphasize the exponential character of 
the PDF tails. 
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Figure 6: PDF of the instantaneous, local fluctuations of the spanwise component of wall 
friction around the time- and space-mean value. The abscissa is in units of standard deviation. 


of such very rare events. The probability is extremely low in all cases, as we are 
speaking of events with a probability of occurrence of 10“® — 10“®. However, 
the probability of observing strong events with cr > 12 is 2-3 times lower for 
CFR than for the other two forcing strategies. 

As shown in Table [U the actual friction Reynolds numbers of the three 
cases are slightly different. In particular, the Reynolds number of CFR is the 
smallest among the three cases, although by a tiny amount. Therefore, one 
might suspect that the present deviation is caused by the small difference in 
Rcr- Indeed, it is known till that there is a slow trend of increase for such 
quantities with Re. However, if this were the case a similar difference would also 
be expected to show up for the the spanwise wall shear stress. Since the PDFs 
of the spanwise wall shear stress as shown in figure |6] do not show discernible 
differences among the three cases, a Reynolds number dependency within the 
results presented in figure [S] may be ruled out. The PDF for the spanwise 
component of the wall friction (figure [6]) are symmetrical due to the reflectional 
symmetry of the problem {z —>■ —z), and the three curves for the different 
forcing strategies basically overlap. Considering that the forcing term in CFR 
shows the largest fluctuations as shown in figure [U it might be expected that 
this directly increases the fluctuations of the spatially-averaged wall shear stress. 
However, figure [5] indicates that the extreme events of the very lare wall shear 
stress are less probable in CFR. Hence, the difference cannot be explained as 
a direct consequence of the forcing and therefore suggests the existence of an 
indirect dynamical effect on near-wall turbulence. 
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Figure 7: Space-time autocorrelation function R of the wall friction at zero spanwise separa¬ 
tion. Lines show levels of the autocorrelation functions from 0.9 down to 0.1 by 0.1 decrements. 
The color code is the same as used before; i.e. red=CFR, black=CPG, green=CPI. 


We now move to two-point, two-times statistics, and show in figure [3 the 
autocorrelation R of the wall friction as a function of the spatial (longitudinal) 
and temporal correlation, computed at zero spanwise separation. The definition 
of R is 


K(Ax.Al) = + + 

'^wrms 


(14) 


where a prime denotes a fluctuating component from the time-space average. 
The autocorrelation functions show the well-known ’’banana” or ’’cigar” shape, 
that indicates the convective nature of the flow: the slope of the inclined narrow 
region, which dimension-wise represents a velocity in such a plane, broadly 
corresponds to the convection velocity of the friction fluctuations [sl ll2| whose 
value is around 10 when expressed in inner units. The autocorrelation of the 
wall friciton has been looked at in the past in several papers, but mainly as a 
function of time separation alone, or space separation alone, i.e. along one of 
the axes of figure [7l Indeed, temporal correlations are the obvious choice for 
experimentalists; spatial correlations are easily computed from DNS, whereas 
exploring the correlation in the space-time plane requires the database to be 
fully time resolved. The three correlation functions computed here are almost 
identical in the whole small-separation region where R assumes large values, say 
R > 0.2. However, the contours at i? = 0.1 do indeed show some difference that 
could be ascribed to something more than residual error. 

This difference was much larger when only a portion of the database was 
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Figure 8: Space-time autocorrelation function R of the wall friction at zero spanwise separa¬ 
tion. Thick continuous lines are computed by using the full statistical sample, dashed lines are 
computed by using the first one-third of the dataset, and thin continuous lines are computed 
using the first two-thirds of the dataset. The dot indicates where the convergence process of 
the autocorrelation function is monitored (see next figure [9jl. The color code is the same as 
used before; i.e. red=CFR, black=CPG, green=CPI. 



Figure 9: Convergence history of the correlation function R for Ax = 8.6 and At"*" = 150. 
The color code is the same as used before; i.e. red=CFR, black=CPG, green=CPI. 
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Figure 10: Two-dimensional autocorrelation of the spanwise component of wall friction (left) 
and pressure (right) at zero spanwise separation. The color code is the same as used before; 
i.e. red=CFR, black=CPG, green=CPI. 

available. To investigate this, a simple visual estimate of how the curves in 
the high-separation region are converging upon increasing time average can be 
obtained in hgurelH where the same correlation function of figure [3 is replotted 
by only using a partial chunk of the whole database. It becomes evident that 
the contours at i? = 0.1 present some residual variations upon increasing the 
statistical sample, and the possibility exists that they eventually converge upon 
further increasing the size of the sample. This view is better supported by 
figure m where the convergence process of the value of the autocorrelation for a 
fixed separation of At+ = 150 and Ax = 8.6 (shown with a black dot in figure 
B is monitored as the computing time progresses. The very slow convergence 
process at such large temporal separations can be appreciated, together with 
the tendency of the three curves to progressively approach the same value. 

To conclude, we show in figure [10] the same spatio-temporal correlations for 
the other quantities that are defined at a solid wall, i.e. the spanwise component 
of the wall shear (shown in the left plot) and the pressure (right plot). The 
contours for both quantities are qualitatively similar to those for the streamwise 
component of wall shear, with the elongation implied by the convective nature 
of the flow and the convective velocity scale shown by the mean slope of the 
contours. The typical Lagrangian timescale associated to these correlations is 
shorter than that for the longitudinal friction, and the convective velocity scale 
is larger, especially for the pressure correlation. What is most relevant here, 
however, is that neither function shows differences that can be attributed to the 
particular forcing strategy. Indeed, the spanwise friction shows some differences 
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at the lowest correlation level R = 0.1, that may be due to a still incomplete 
convergence as in the case of longitudinal friction. 

4. Concluding remarks 

Since the choice of a specific forcing term in the DNS of a turbulent channel 
flow is to some extent arbitrary, in this paper we have searched for possible 
statistical consequences of this choice. It is not inconceivable that such signature 
exists, as different forcing terms produce different dynamical systems, unless an 
infinitely large computational domain is considered. 

The difference between the forcing terms for CFR, CPG and CPI can indeed 
be seen in the time traces of instantaneous global flow properties. In order 
to obtain an indication whether these differences are also reflected in some 
statistical properties of the flow, three databases have been built by DNS, one 
for each forcing strategy. 

As far as one-point one-time statistics are concerned, the agreement among 
the results of the three forcing strategies is extremely good. This was the ex¬ 
pected outcome, as extensive DNS research over several decades did never iden¬ 
tify any effect of the forcing term in such frequently encountered statistics. 
However, we have observed small yet statistically significant differences in the 
PDF of the streamwise wall friction, where CFR shows fewer very rare events 
of extremely high wall friction. For such differences to become evident, a very 
large database is required as the PDF begin to deviate for very rare events, 
whose intensity is cr > 12. In other words, in such events the amplitude of the 
friction fluctuations above the mean is almost 5 times the mean value, and their 
probability of occurrence is about I0“®. Small differences might also be hidden 
in the space-time autocorrelation of wall friction, although this quantity is more 
sensitive to statistical error and the size of our database is still not large enough 
to enable a firm conclusion on this matter. 

Hence, we have found no reason for questioning the way DNS research of 
turbulent flows is carried out, i.e. by choosing whatever forcing term turns 
out to be the most convenient in practice. However, having found at least 
one statistical quantity (the wall friction PDF) that is affected by the specific 
forcing term is an important result in principle. It can also be envisaged that 
other quantities too might carry traces of the forcing term, or that differences 
might be observed in a Lagrangian framework. In general, these differences 
are likely to be more pronounced as the size of the computational domain is 
reduced. For sure, they matter whenever a flow is changed from one state to 
another, which is typical, for example in flow control. The choice of the forcing 
term might also be important in studies on laminar to turbulent transition. 
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